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ABSTRACT 

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of 
its shape if the PSF profile depends on wavelength. We study how this bias depends on the 
properties of the PSF and the galaxies themselves. The bias depends on the scales used to 
estimate the shape, which may be used to optimise methods to reduce the bias. Here we 
develop a general approach to quantify the bias. Although applicable to any weak lensing 
survey, we focus on the implications for the ESA Euclid mission. 

Based on our study of synthetic galaxies we find that the bias is a few times 10 -3 for 
a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be 
accounted for. We demonstrate how one can do so using spatially resolved observations of 
galaxies in two filters. We show that HST observations in the F606W and F814W filters al- 
low us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's 
scientific requirements. The precision of the correction is ultimately determined by the num- 
ber of galaxies for which spatially-resolved observations in at least two filters are available. 
We use results from the Millennium Simulation to demonstrate that archival HST data will be 
sufficient for the tomographic cosmic shear analysis with the Euclid dataset. 

Key words: Gravitational lensing:weak, surveys 



1 INTRODUCTION 

The measurement of the distortion of the shapes of galaxies caused 
by gravitational lensing by large scale structures, i.e. cosmic shear, 
is a powerful tool to investigate the statistical properties of the Uni- 
verse and in particular to understand the mechanism responsible 
for the observed accelerated expansion. The origin of this accelera- 
tion has been dubbed "dark energy" reflecting the fact that we still 
lack a theoretical framework to explain its nature. To discriminate 
between competing theories, large cosmological experiments are 
needed. The most ambitious among these is the recently selected 



ESA mission, Ewc/i'4]] fcaureiis et alj|201lh . which will survey the 
15,000 deg 2 of the extragalactic sky that has both low extinction 
and zodiacal light. 

The measurement of the second-order shear statistics as a 
function of redshift, which is commonly referred to as cosmic shear 
tomography, is particularly powerful to study th e growth of st ruc- 
tures and the expansion history of the Univ erse dHdl999ll2002h see 
Hoeks tra & Jainl2008l ; lMunshi et alj2008l. for recent reviews). The 
accuracy of the derived constraints on the cosmological parameters 
depends critically on our ability to measure both the shapes and 
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redshifts of the billions of galaxies that will be observed. Both are 
challenging tasks and Euclid has been designed to achieve this. For 
a detailed discussion o f how this can be done, we refer the reader 
to lCropper et alj d2012h . 

A critical step in the estimation of the cosmic shear sig- 
nal from the observed shapes of galaxies is the correct i on for 
the point spread function (PSF): as shown in lMassev et al.1 d2012h 
most residual biases scale proportional to the square of the 
PSF size. The biases affecting shear measurement techniques 
have been studied in a number of colla borative projects such 
as the Shear TEst ing Programme (STEP; iHevmans et ail 120061 : 
iMassev et alj|2007|). the Gravitational LE nsing Accuracy Testing 
2008 ( GREAT08;lBridle et alj20"09L|20ld) and the GREAT 10 chal- 
lenges dKitching et alj|20ict l2012a lbl) using simulated monochro- 
matic data. The impact of the noise in the data, which affects the 
fainter galaxies, has been recently studi ed (Refregie r et alj|2012l ; 
iKacprzak et"al]|2012l ; iMelchior & Viobj|2012l) . While these stud- 
ies are necessary to calibrate methods so that they will perform in- 
creasingly well on monochromatic images, we explore here another 
source of bias which has not been considered by those studies. 

Cosmic shear surveys are designed to maximise the number 
of observed galaxies. For this reason they use broad-band filters; 
the observed images are the integrated light distribution over a 
large wavelength range. The spectral energy distribution (SED) of 
a galaxy typically varies spatially, generating "colour gradients". 
This prevents one from unambiguously recovering the unconvolved 
light distribution required for an unbiased shear estimate from the 
observed images. 

The existence of colour gradients is one potential source of 
systematic error diminishing the power of future weak lensing mis- 
sions. Observations of galaxies in a broad-band filter to estimate 
the shear add some complications even if galaxies have no colour 
gradients. This happens because the PSF is chromatic. Ignoring thi s 
fact leads to a bias which was quantified bv lCypriano et alj J2010h . 
This bias can be avoided using a colour weighted PSF, that depends 
on the global SED of the galaxy. In the case of Euclid this PSF can 
be estimated with sufficient accuracy from supporting deep multi- 
colour ground-based observations. 

While the PSF correction described bv lCypriano et alj d2010l) 
is perfect when the galaxies have no colour gradients, it becomes 
inaccurate when the colour varies spatially. This can be understood 
by considering the simple case of a galaxy with a small red bulge 
and a much more extended blue disk. The global colour that is used 
to estimate the SED-weighted PSF will resemble that of the disk if 
its flux exceeds that of the bulge. To suppress noise, shape measure- 
ments employ a radial weighting that is matched to the brightness 
profile of the galaxy. This enhances the contribution of the bulge to 
the shape measurement, implying that a PSF redder than the SED- 
weighted one should be used. 

As we will quantify in detail below, the bias depends on the 
width of the filter that is used. Consequently, colour gradients are 
expected to be particularly rele vant for Euclid because of its wide 
pass-band dLaureiis et alj|201 lh . We note, however, that it may not 
even be negligible for future ground-based experiments. Despite 
the fact that multi-wavelength ground-based projects will employ 
narrower filters and that the wavelength dependence of the PSF 
is weaker, the intrinsi cally larger PSF e xacerbates the impact of 
colour gradients (see lMassev et al.ll2012t for a detailed analysis of 
the impact of the PSF size and other systematics on shape measure- 
ments). 

The impa ct of colour gradi ents on shear measurements was 
first studied bv lVoigtetal] j2012) who used a small sample of HST 



alaxies for which a bulge plus disk d ecomposition was available 
Simard ct al]2002h . lVoigt et alj J2012h quantified the bias induced 
by colour gradients and also estimated the minimum size of a sam- 
ple of galaxies that needs to be observed in two narrower filters in 
order to determine the bias. Because of the small sample size and 
a concern that the bulge plus disk de composition mig ht not fully 
capture the properties of real galaxies, IVoigt et alj 1 2012) provided 
conservative upper limits. 

The aim of this paper is to expand upon this pioneering study. 
To this end we develop an approach that allows us to explore how 
the bias depends on the wavelength range covered by the filter and 
on the properties of the PSF and the galaxies. Furthermore, we dis- 
cuss how the amplitude of the bias depends on the method used to 
estimate the shear. To do so, we use simulated galaxies and show 
how the bias can be modelled. Our methodology can in principle be 
applied to archival HST data or any other dataset of well-resolved 
galaxies observed in least two filters. We argue that the HST archive 
represents the most suitable dataset for colour gradient studies be- 
cause of its PSF characteristics. We show that a sufficient number 
of galaxies has been observed to calibrate the colour gradient in- 
duced bias with the precision required to achieve Euclid's science 
objectives. 

The paper is organised as follows: in Section[2]we introduce 
the notation, we describe the nature of the problem and the way we 
quantify the bias. In Section[3]we describe how we produce simula- 
tions which we then use in Section|4]to study the size of the bias as 
a function of the characteristics of the PSF and of galaxies. In sec- 
tion|5]we show that it is possible to construct a calibration sample 
to model the bias using observations with two narrower filters. In 
Section[6] we discuss the performance one can achieve using obser- 
vations in the F606W and F814W HST filters. We show in Section 
[7] that the number of galaxies observed in the F606W and F814W 
(or F850LP) in the HST archive will be large enough to characterise 
the bias with the required precision. We conclude in Section[8] 



2 DESCRIPTION OF THE PROBLEM 

We start by showing how the spatial variation of the SED across a 
galaxy profile affects the shear estimation and introduce the nota- 
tion that will be used throughout the paper. To facilitate the read- 
ability of the paper we summarise the main quantities in Table[T] 

Throughout the paper we implicitly assume that measure- 
ments are done on images produced by a photon counting device, 
such as a charge-coupled device (CCD). Hence J(x;A), the ob- 
served photon surface brightnesjf| or image at a wavelength A, is re- 
lated to the intensity 5(x; A) by /(x; A) = AS(x; A)T(A), where 
T(A) is the normalised transmission. The image of a galaxy ob- 
served with a filter of width AA is given by 



I obs (x)= / /°(x;A)*P(x; 



\)d\ . 



(1) 



where P(x; A) is the wavelength dependent PSF, and /°(x, A) is 
the image of the source before the convolution (denoted by *) with 
the PSF. In Fourier space, Equation. (fTJ can be written in a more 
convenient way as 



2 Note that we drop the explicit mention of 'photon' after the first part of 
Section 2. 
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/ obs (k) = / 7°(k;A)P(k;A)dA. 



(2) 



To measure cosmic shear one needs to estimate the second-order 
moments of the PSF-corrected image 7° ( x) = f A a 7°(x; A)rfA. In 
the w eak lensing regime (see for example Bartelmann & SchneideJ 
l200lh the shear can be estimated from the measurement of second- 
order moments Qy : 



7i + 172 



Q11 - Q22 + HQ 



Q11 + Q22 + 2(Q f 1 , 1 Q22 



,2-11/2 



(3) 



where we have defined the complex shear 7 = 71 + 172. The 
second-order moments of the light distribution are given by 



Q" 



dX I I (x; X)xiXjdit , 



(4) 



where we implicitly assumed that they are evaluated around the po- 
sition where the dipole moments vanish. The total observed photon 
flux is given by 



1 f 1 2 

b Ja\ 2 

+~ f (/ifei + \q^a) (A - A ) 2 dA + 0(AA 4 ) , 



' AA 

where we have defined 



fk 



d k F(X) I 



and 



d\ k U ref ' qij ' k ~ dx k U ref ' 

We have split the right-hand side of Equation <[Sj> into two terms on 
purpose. The first term corresponds to the second-order moments if 
the galaxy had no colour gradients. The second term represents the 
lowest order correction due to colour gradients. By evaluating the 
integral we obtain: 



Q n = Q 



nograd 



(AX) 2 f q^h , 1_ 
12 \F{X tei ) + 2 qii ' 2 



0(AX 4 



(9) 



which shows that the change in the quadrupole moments is propor- 
tional to (AA) 2 . Note that 



F = 



f F(\)dX = f dX f J°(x;A)dx. (5) 

J AA J AX J 



The observed moments are measured from the PSF-convolved im- 
age given by Equation (TJ. In practice, to reduce the effect of 
noise in observed images, moments are evaluated using a weight 
function Wdxl) with a characteristic size r w . Hence the observed 
quadrupole moments are given by 



?obs 



1 
"75" 

1 w 



d\ J dx7°(x;A) *P(x;X)xiX j W{\x\: 



(6) 

with F w indicating the weighted flux (i.e. the weight function 
is also introduced in Equation ((5}), from which it is still pos si- 
ble to retrieve the s hear (see for example iKaiser et al.l 1 19951 or 
lMelchioretall201lh . 

The alternative to mome nt-based techniques is to use so- 
called fitting techniques (e.g. iBridle et al.l [20021; iKuiikenl 1 20061 : 
iMiller et al.ll2007l : iKitching et alj|2008l : iMiller et alj|20l3> 7 which 
fit the sheared, PSF-convolved galaxy profiles and thus provide an 
estimator of the shear. In this case, the profile itself acts as a weight 
and determines the parts of the galaxy profile that are used to esti- 
mate the shear. The weighting is a very important point to keep in 
mind because, as we will show, the amplitude of the colour gradi- 
ents induced bias strongly depends on the characteristic scale of the 
weight function. To highlight the importance of the weight func- 
tion, we first examine the impact of colour gradients on unweighted 
moments. We start by noting that the expression for the quadrupole 
moments can be written as 



h 



1 



dF(X) I 



F(A r . 



F(X t 



OX 



dlnF(A) 
dX : 



(10) 



which means that the change of the moments depends on the vari- 
ation of the SED of the galaxy across the filter, i.e. on its colour. 

Interestingly, it is still possible to determine from the ob- 
served unweighted quadrupole moments. This can be seen by writ- 
ing down the convolution explicitly an d changing the or der of inte- 
gration. Doing so we obtain (also see lValdes et al J 19831) : 



<2°j bs = -^ J -?°(x; A) * P(x; X)xiXjdX 



■ -p I F(X)P tJ (X)dX. 



(11) 



The observed quadrupole moments can thus be written as the sum 
of the true moments and Pij(X), the quadrupole moments of the 
PSF integrated over the pass-band. More generally, even moments 
of order N of 7°(x) are expressions involving the moments of 
7°(x; A) up to TV — 2. In the case of the second-order moments, 
only the knowledge of the SED is needed to retrieve Qy. Hence 
colour gradients do not bias the shear estimate based on unweighted 
moments, provided the SED of the galaxy and the PSF moments as 
function of the wavelength, Pij (A), are known. 

In contrast, it is not possible to correct for the PSF without 
knowledge of the higher order moments when weighted quadrupole 
moments are used. It is therefore not possible to recover 7°(x), ex- 
cept in two special cases. The first case is when the PSF is achro- 
matic: 



- i :i 



F 



)ij(X)F(X)dX , 



(7) 



where (A) is given by Equation 01 with an infinitesimally nar- 
row filter centered on A. For the following it is convenient to define 
A rc f, a reference wavelength such that F = F(X re i)AX. Using a 
second order Taylor expansion for both F(X) and Qlj{X) around 
A rc f , and keeping only the even powers (the odd powers will vanish 
after integration) one finds that 



The second case is when the galaxy has no colour gradients: 



7°(k; A rof ) 



F(A rcf )7 obs (k) 
J AA F(A)P(k;A)dA' 



(12) 



(13) 



and one can rewrite 7(x; A) = 7°(k; A rG f )F(X) / F(X Te t) for any 
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Symbol Description 





Second order moments before PSF convolution. 

Without explicit A dependence it refers to the moments integrated over the whole pass-band. 


7°(x) 


Photon surface brightness or image describing the source light distribution before smearing by the PSF. 
Without explicit A dependence it refers to the counts integrated over the whole pass-band. 


S(x; A) 


Source intensity as a function of position and wavelength. 


Q° hs 


Observed second order moments after PSF convolution. 

Without explicit A dependence it refers to the moments integrated over the whole pass-band. 


7" obs (x) 


Observed photon surface brightness or image. 

Without explicit A dependence it refers to the counts integrated over the whole pass-band. 


PijW 


Second order moments of the PSF. 

Without explicit A dependence it refers to the moments integrated over the whole pass-band. 


Rpsf 


Characteristic size of the PSF estimated as the sum 
of the second order moments: V-Pti + P22 




Characteristic size of the observed galaxy estimated as the sum 
of the second order moments: -J QJ bs + Q^ s 



T(A) Transmission as a function of wavelength A. 



'y Estimate of the complex shear vector. 



Qnograd Second order unconvolved moments of a galaxy with no colour gradients. 



F(X) Photon flux at a given wavelength A, such that F = f F(X)dX = f \S(X)dX, where S is the flux density. 



Observed complex ellipticity. 



Response of a galaxy ellipticity to a shear 7. P" ogl ad is the response of a galaxy that is assumed to have no colour gradients, 



while P grad is the true response. 



/J o rad (x) Image constructed assuming a galaxy has no colour gradients. 

This quantity is needed to estimate the bias as described in Section lTTI 



^nograd( x ) Observed image obtained after applying a shear to 7^ ograd (x) and convolving by the PSF. 

Without explicit A dependence it refers to the moments integrated over the whole pass-band. 

7 ob ad (x) Observed image obtained after applying a shear to the true galaxy profile / (x; A) and convolution by the PSF. 
Without explicit A dependence it refers to the moments integrated over the whole pass-band. 

Table 1. Summary table of quantities defined in this paper. 



given choice of A rc f. In this last case one is able to derive the PSF 
corrected image by knowing both the flux of the galaxy and the 
PSF profile as a funct ion of wavelength. This is the case studied by 
ICvprianoet all ( bOld) . 



2.1 Measurement 

The need to use a weight function when measuring the shapes of 
the faint galaxies leads to a bias in shear estimates due to colour 
gradie nts. This result ap plies to all mo ment-based methods, such as 
KSB dKaiseretalJl995t) or DEIMOS jMelchior et alj201 ll) . How- 
ever, the fact that one cannot recover the unconvolved image 7° (x) 
from broad-band observations suggests that 'fitting methods' will 
also be prone to provide biased estimators of the shear unless they 
are able to account for the existence of colour gradients. For in- 



stance, one could attempt to model galaxies with a bulge and disk 
component, each with their own SED. We proceed to quantify this 
bias for moment-based methods, but note that our approach can be 
extended to evaluate the bias for fitting techniques as well. 

The first step is to define the bias that is induced by a spatially 
varying SED. The PSF correction described by Equation J 1 3 b is 
perfect for a galaxy that has no colour gradients. If Equation jl3\ is 
used to obtain the unconvolved image of a galaxy, one effectively 
approximates the observed galaxy with a galaxy that has the same 
SED, but no colour gradients. This galaxy actually has a different 
profile and thus its response to the shear and to the PSF is different. 
For this reason, correcting the galaxy using Equation J 1 3 fc leads to 
a biased estimate of the shear. To quantify the resulting bias we 
define the response P 7 as the link between the observed ellipticity 
e obs and the shear 7 
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Figure 1. Spectral energy distributions used to create the disk and bulge 
components of our synthetic galaxies. All SEDs are normalised such that 
the integrated S"(A) between 200 and 1100 nm is one. In our synthetic 
galaxies the SED of the bulge (red solid line) is fixed, while the SED of the 
disk can be either Sa, Sb, Sc, Sd simulating a bluer and bluer galaxy, or Irr 
which represents a starburst population. All reference galaxies (see Table[5J 
use an Irr SED. 



(14) 



where we defined the observed complex ellipticity using weighted 
quadrupole moments: 



+ ie 2 = 



-jobs 

hi 



Qobs | r>-/-)obs 
22 + 2tQl2 



Qobs + Qobs + 2 (Q ? bsQ° 



(0^)2)1/2 



(15) 



In the absence of observational biases, Equation d 1 4b can be used 
to obtain an unbiased estimate of the shear. In practice the shear is 
obtained by averaging over many galaxies, since they have an in- 
trinsic ellipticity that is much larger than a typical shear. In the case 
of a galaxy with colour gradients one is unable to estimate P grad , 
the correction to apply to obtain an unbiased shear estimate. One 
will instead correct the PSF using Equation l |13| > which replaces the 
response P grad with P" ogra . This approximation leads to a mul- 
tiplicative bias which we define as: 



21-i = 

7< 



^>nograd 



(16) 



where ji refers either to the first or second component of the shear 
pseudo-vector; we make the assumption that the bias m is the same 
for both components. To estimate the difference in response we cre- 
ate a pair of galaxies which appear identical when observed through 
a broad pass-band. One of them has a colour gradient and the other 
does not. We apply the same shear to both galaxies, convolve them 
with their respective PSFs and derive pn°s rad an d P| rad by com- 
paring the observed ellipticities to the applied shear. 



to determine the bias using real data. Here we examine instead 
whether it is possible in principle to measure the colour gradient 
bias with sufficient precision for Euclid. We do so by making a 
number of conservative assumptions, while using synthetic galax- 
ies and an analytic description of the PSF. 

We assume that galaxies can be described as the sum of a 
bulge and a disk component, each characterised by Sersic profiles 
of index n. The Sersic profile of index n is given by: 



S(x) = S c e 



- fC [(x-X ) T C(x-X )]2r l 



(17) 



where x T is the transpose of x, S c is the value of the intensity in 
the ce nter xq, k = 1.9992n — 0.3271 (see for example [Capacciolil 
1989), and C is a matrix defined by: 



f cos 2 Q) sin 2 (0) 

2/JA 2 



C22 = (™-. 

\ a 



sin ((f)) cos (4>) 
T~ + 7j2 



(18) 
(19) 
(20) 



with <f> the angle between the semi-major axis a and the x-axis and 
b the semi-minor axis. In the case of an axisymmetric profile ru = 
a = b is the half-light radius. The ellipticity of the profile is defined 
as e = (a — b)/(a + b). 

We acknowledge that this double Sersic model may not de- 
scribe the full variety of observed galaxies. As our aim is not to 
determine the bias, but rather develop an approach to measure the 
bias, this model is sufficient for our purposes as it should cap- 
ture most of the problem of colour gradients. Whenever statements 
which depend on the value of the bias have to be made, we will 
make sure that our assumptions are conservative and the bias is not 
underestimated. 

To describe the S ED of the bulge and the di s k we use the 
galaxy templates from IColeman. Wu. & Weedmanl (l980) shown 
in Figure [TJ The bulge is always modelled with the SED of an old 
stellar population typical of elliptical galaxies, while the disk is 
modelled either with an extremely blue SED typical of an irregular 
galaxy, or intermediate stellar populations such as Sa, Sb, Sc and 
Sd. For each wavelength we construct the galaxy surface brightness 
profile 5°(x; A) by adding up the profiles of the bulge and the disk 
normalised by fixing the ratio Sbuigc/5'disk at A = 550 nm. The 
profile is sampled in a wavelength range between 200 and 1100 
nm with a step of 1 nm. The pixel size of each S° (x; A) image is 
0.05 arcsec. 

In order to generate the observed synthetic galaxy image we 
need to decide on the shape of the transmission curve T(A) and 
the profile of the PSF at any given wavelength. The black solid 
line in the left panel of Figure [2] shows the curve we use to model 
the expected transmission for Euclid. Note t hat we define here a s 
EuclidPSF, the PSF used in the Design Study dLaureiis et al.l20lTh . 
This PSF is not the final one but we do not expect major changes. 
We use as reference a PSF with a Gaussian profile and wavelength 
dependent width: 



3 SIMULATIONS 

The main aim of this paper is to develop a method that can be used 
to estimate the bias induced by the presence of colour gradients. In 
a future paper (Huang et al. in prep.) we will apply our approach 



a(A) = ^ 800 (8TOn^) a ' (21) 

where a is the sl ope and wp.soo is the width of t he P SF at 800 nm. 
As discussed bv lPaulin-Henriksson et alj d2009l) and lMassev et alj 
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Figure 2. Left panel: The adopted Euclid transmission curve T(A) (black solid line) and those for the HST filters F606W (blue solid line) and F814W (red 
solid line) used in Section[3]and[6]to model the bias. Note that the normalisation of the Euclid transmission curve is arbitrary, whereas for the F606W and 
F814W filters we show the measured throughput T(X). Right panel: comparison of the Euclid PSFs at 550 nm (solid blue line) and 800 nm (solid red line) 
and our reference model PSF (PSF1; dashed blue and red lines). The parameters for our model PSFs are listed in Table|2] PSF1 is broader and has a stronger 
wavelength dependence than the actual Euclid PSF. The dotted lines show the profiles for PSF3, which resembles the Euclid PSF more closely (see also 
Sectionl4~2l. 



d2012l) . the ability to recover the shear of a galaxy of given char- 
acteristic size 7?g al = Q°i s + Q§2 S is oc -Rp SF /-Rg a i where 
i?P SF = Pn + P22 is the characteristic size of the PSF. This re- 
flects the fact that the PSF blurs the galaxy and this effect can be 
only partially corrected. 

We construct our reference PSF by choosing a = 1 and 
wo, 800 = 0.102 arcsec. This yields a PSF with the same char- 
acteristic size -Rpsf as the Euclid PSF for A = 800 nm but with a 
stronger wavelength dependence (see right panel of Figure^. For 
comparison purposes we also define an even broader PSF, indicated 
as PSF2 in Table [2] The Euclid PSF is not well approximated by a 
single Gaussian, but instead can be considered as being composed 
of two parts: the core with a diffraction limited characteristic size 
oc A and the wings which also contribute to Py but have a much 
weaker wavelength dependence. For this reason the overall Euclid 
PSF size scales as P PSF oc A ' 55 . Its FWHM does scale propor- 
tionally to A as it depends on the core. 

Our reference PSF is much broader than the Euclid PSF and 
this leads to bias estimates considerably larger than the ones ex- 
pected for the actual Euclid PSF. To quantify how conservative 
the results are when using the reference PSF, we compare to PSF3 
which is the sum of a Gaussian and a top-hat. The Gaussian, with a 
width luo.soo = 0.054 arcsec, is chosen to fit the core of the Euclid 
PSF, while the top-hat approximates the wings. The cut-off size of 
the top-hat is approximately at the position of the second minimum 
of the PSF profile, r cut _ ff oc A 0,74 and the normalisation is such 
that the top-hat contains 20% of the total flux. As shown in the 
right panel Figure[2] PSF3 approximates the main properties of the 
Euclid PSF fairly well. 

Using the model for the galaxy, the PSF and the transmission, 
we compute the observed image J° bs (x). We proceed following the 
procedure outlined in Section|2]and create pairs of galaxies that ap- 
pear identical in the broad-band image, but where one has a colour 
gradient and the other does not. To create the latter, we deconvolve 
7 obs (x) using the colour- weighted PSF correction in Equation H3\ 
and call this /no g rad( x i We apply the same shear to J°(x; A) 
and /nograd( x ; -*0- We convolve the profiles by the PSF P(x; A) 



PSF Description 

Reference (PSF1) Gaussian PSF described by Eqn. )2U 

with too,800 = 0.102 arcsec and a = 1. 

Wide Gaussian (PSF2) Gaussian PSF described by Eqn. <2 1 1 . 

with i«o,800 = 0.15 arcsec and a = 1. 

Gaussian + top-hat (PSF3) Gaussian core described by Eqn. i2\\ 

with i«o,800 = 0.054 arcsec and a = 1 
and top-hat with 20% of the total flux 
and a cut-off size ex A 0,74 . 

Table 2. Summary table of the PSFs considered in Section 14.21 The refer- 
ence PSF is used throughout the paper, as its large width exacerbates the 
bias. PSF3 resembles the actual Euclid PSF more closely. 

and create the final observed images by summing /g r ad(x; ^) an d 
J°og rad (x; A) over the full pass-band. Note that we use a pixel size 
of 0.05 arcsec, because we do not want undersampling to affect the 
estimate of the bias. To construct the Inograd an d -fgrak images one 
would want the best resolution available to avoid sampling bias. 
The sampling we use corresponds to the HST ACS/WFC resolu- 
tion, which is sufficient to avoid undersampling of the Euclid PSF. 
We then measure the ellipticities from these images using weighted 
quadrupole moments and determine P 7 for each galaxy. To reduce 
noise in our estimate of the multiplicativ e bias we use the 'ring- 
test' method (Nakaiim a & Bernsteinl2007l) creating 8 copies of the 
same galaxy but with different orientation. 



4 EVALUATION OF THE BIAS WITH SYNTHETIC 
GALAXIES 

The intrinsic properties of the source galaxies vary with redshift, 
and the colour gradient bias will therefore differ b etween tomo- 
graph ic redshift bins. As shown below (but also see lMassev et alj 
|2012|) the bias depends on the intrinsic galaxy size and for this rea- 
son we define a number of model galaxies in Table [3] that span a 
range in size. The properties of these galaxies are chosen purposely 
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Name 


SED 


a (arcsec) 


Flux ratio (550 nm) 


Sersic index 


FWHM (z=0)(arcsec) 


FWHM(z=0.9) 


B 


E/Irr 


0.17/1.2 


25%/75% 


1.5/1.0 


0.32 


0.40 


S 


E/Irr 


0.09/0.6 


25%/75% 


1.5/1.0 


0.27 


0.37 


B4 


E/IlT 


0.39/0.35 


33%/67% 


4.0/1.0 


0.22 


0.30 


S4 


E/Irr 


0.12/0.20 


33%/67% 


4.0/1.0 


0.22 


0.27 



Table 3. Characteristics of our reference galaxies. When two values are quoted in a column, the first value refers to the bulge the second to the disk. The 
galaxies are circular (a = b) with the size of the semi-major axis a indicated in the second column. The observed FWHM at z = and z = 0.9 are 
computed using the reference Gaussian PSF. Note that for these galaxies we keep the Sersic parameters of disk and bulge fixed as a function of redshift. As a 
result, the FWHM increases at high redshifts because the dis k becomes b righter tha n the bulge. A secon d pair of galaxies (B4 and S4) is chosen such that the 
bulge-to-disk properties are similar to the typical galaxy from lVoigt et"ai] <2012h and lSimard etsdl <2002h . 



0.0100 r\ 



0.0010 



0.0001 




Figure 3. Amplitude of the absolute value of the bias as a function of r w , 
the width of the weight function, in units of r h , the half-light radius, for the 
reference galaxies B (solid black line) and S (dashed black line) defined in 
Tablefjjfor 2 = 0. The dotted red line indicates where the weight function 
size is equal to rh . The blue solid (dashed) line indicates the noise-to-signal 
ratio for the B (S) galaxy as a function of the filter size normalised to the 
value obtained when r w = r^. Note that the noise-to-signal ratio shown 
here only includes the contribution from the sky background (see text). 




— 0.004 



0.000 



0.0 0.2 



Figure 4. The top panel shows the value of the bias for the B galaxy as 
a function of redshift for the PSFs in Table [2] The black solid line shows 
the result for the reference PSF (PSF1); the red solid line shows the bias 
for a Gaussian with a larger width (PSF2), while the blue line indicates the 
bias when the PSF is approximated by a Gaussian with width mio,800 = 
0.05 arcsec plus a top-hat (PSF3) whose profile is shown in right panel of 
FigurefS] The bottom panel shows the same results but for the S galaxy. 



to have large colour gradients: they are a superposition of a rather 
bright red bulge with small characteristic size and an extended disk 
with an SED of an irregular galaxy. 

We note that the reference PSF used to compute the observed 
FWHM values listed in Table [3] is much broader than th e Euclid 
PSF (see Fig. [2}. As discussed in lLaureiis et alj d201lh the Eu- 
clid source galaxy sample is selected to have an observed FWHM 
> 1.25 x FWHMpsf (where FWHM PS f ~ 0.15" is the FWHM 
of the PSF at A = 800 nm). Comparison to observed sizes in deep 
HST data indicate that galaxy S is actually representative for the 
smal lest galaxies that w ill be used in the Euclid weak lensing anal- 
ysis JMassev et alj|2012l) . 

4.1 Dependence on the weight function 

As mentioned in Section [2] the bias is a function of the weight 
function, which we take to be a Gaussian with a dispersion r w . The 
optimal choice, in terms of maximising the signal-to-noise ratio, 
is to match the weight function to the size of the source galaxy. 
Figure [3] shows the amplitude of the bias as a function of r w for 
the B and S model galaxies defined in Table |3] For the reference 
galaxy B we obtain a bias of approximately —8 x 10~ 4 if we take 
r w to be equal to the half-light radius rh, which was measured on 
the original observed image /g£fd( x ) (i.e., without applying shear). 



Compared to the B galaxy, which has the same bulge-to-disk flux 
ratio, the colour gradient-induced bias is larger for the S galaxy 
because its bulge and disk sizes are both a factor of two smaller. In 
this case the bias is —3 x 10~ 3 when r w = rh- 

The bias increases when r w decreases, whereas the bias van- 
ishes when r w goes to infinity, as expected from SectionfJ] This re- 
sult can be understood by noting that the PSF correction described 
by Equation J 1 3 b is weighted by the colour of the galaxy. For galaxy 
B, the disk contains 75% of the flux whereas the bulge contains 
only 25% of the total flux. When r w is small, the moments are 
measured essentially from the profile of the bulge, and corrected 
using an effective colour of the PSF which is always wrong. For 
the disk the colour-weighted PSF is closer to the correct one and 
therefore the bias is reduced when r w increases. 

Throughout the rest of the paper we will choose the value 
for r w to be equal to the observed galaxy half-light radius rn be- 
cause this is the optimal choice in terms of signal-to-noise ratio, 
and therefore routinely used in practice. This choice, however, may 
no longer be optimal when colour gradient-induced biases are con- 
sidered. To show this, we compare the value of the bias with the 
noise-to-signal ratio of the observed ellipticity as a function r w . We 
assume that the pixel noise is dominated by the sky background, 
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Figure 5. Amplitude of the absolute value of the multiplicative bias for the B galaxy described in Table [3]as a function of redshift, when varying the bulge 
and disk parameters and using PSF1 (see TablefS). Note that the Sersic parameters of the disk and bulge are fixed as a function of redshift. For comparison, in 
each panel the black solid line indicates the bias for the reference galaxy B. Top panels: from left to right we vary the percentage of the light in the bulge, the 
semi-major axis value of the disk (in arcsec), adisk. the Sersic index of the bulge, the ellipticity of the disk. In the bottom panels we vary: the characteristic 
size of the bulge abulgc ( m arcsec), the Sersic index of the disk, the ellipticity of the bulge, the SED of the disk. 



which is the case for the majority of galaxies in a typical weak 
lensing survey. 

The blue lines in Figure[3]indicate the noise-to-signal ratio in 
the measurement of the ellipticity as a function of r w . Since we 
are only interested in the relative change, the noise-to-signal ra- 
tio is normalised to its value when r w = rn. Because we ignored 
the contribution to t he noise that d epends on the galaxy light profile 
(see Appendix A of lHoekstra et al.l2000f) the ratio is approximately 
constant for small r w . Including this term would result in an upturn 
in the noise-to-signal ratio for small r w , such that taking r w ~ m 
minimises the noise-to-signal ratio. Figure[3]shows that the colour- 
gradient bias decreases rapidly when r w is increased. The noise-to- 
signal ratio increases but does so relatively slowly (note the loga- 
rithmic scale for the bias and the linear scale for the noise-to-signal 
ratio). This suggests that adopting r w > rh might provide a good 
compromise between reducing the amplitude of the bias and in- 
creasing the noise-to-signal when one needs to account for colour 
gradients. 

4.2 Dependence on the PSF and galaxy characteristics 

In this section we study how the bias depends on the PSF charac- 
teristics. The bias as a function of redshift for the two reference 
galaxies B and S is shown in Figure [4] The black solid line indi- 
cates the result obtained using the reference PSF (PSF1 in Table 
O. The red solid line shows the result for PSF2, which has a larger 
width resulting in an increase in the bias. More realistic estimates 
for the actual bias for Euclid are obtained using PSF3, indicated by 
the blue lines in Figure|4j the bias is a factor ~ 2 smaller compared 
to the results for our reference PSF. 

The results presented in Figure [4] are obtained by changing 
the redshift, but keeping the input Sersic parameters of the galaxies 



fixed. This leads to a change in the observed FWHM because the 
disk becomes brighter than the bulge as the redshift increases. Since 
the simulated galaxies have rather extended disks their observed 
sizes also increases (see the comparison of the FWHM at z = 
and z = 0.9 listed in Table[3}. Hence, the value of the bias changes 
with redshift in part because the ratio of the galaxy size to the PSF 
size changes. 

Many parameters contribute to the actual spatial colour varia- 
tions and it is therefore useful to examine how the bias changes as 
a function of the bulge and disk characteristics. To do so, we take 
galaxy B and vary one parameter at a time, keeping all others con- 
stant (using PSF1 and adopting r w — r^). The results are presented 
in Figure[5] We find that the bias depends most strongly on the flux 
of the bulge, but also depends on the size of both bulge and disk. It 
does not depend much on their ellipticity . 

In general the bias values range between values of a few times 
1CP 4 to a few times 10~ 3 . Note that in all panels the rest-frame 
colour of the galaxies is the same, except when we vary the frac- 
tion of the flux in the bulge (top left panel) or when we change 
the SED (bottom right panel). These results confirm what we con- 
cluded based on Equations J9lll0b, i.e. the bias is foremost a func- 
tion of the colour of the galaxy. Whereas the absolute value of the 
bias increases when the size of the source galaxy decreases, the 
changes as a function of the parameters are similar for all galaxies 
in Table[3] The reference galaxy and its smaller version both have a 
relatively large disk, which leads to a large FWHM at z = 0.9 (see 
Table [3}. Therefore we also consider two galaxies with a smaller 
FWHM: B4 and S4. The ratio FWHM/FWHM PSF at z = 0.9 is 
1.3 and 1.1 for B4 and S4, respectively. As one can see from Figure 
[6] the bias is a very strong function of this ratio: it is smaller than 
10~ 3 for a well-resolved galaxy and can become a few percent for 
a galaxy which is about the size of the PSF. 
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Figure 6. Absolute value of the multiplicative bias as a function of the ra- 
tio of the galaxy FWHM to the FWHM of the PSF for the model galax- 
ies S4, B4, S and B (from left to right). The results are for z = 0.9, 
which corresponds to the median redshift of Euclid. The black points in- 
dicate results for the reference Gaussian PSF (PSF1 in Table |5J for which 
PSFpwHM = 0.241 arcsec. Blue values show the results for PSF3 
(PSFpwHM = 0.130 arcsec), which is a better approximation of the Eu- 
clid PSF. Dotted lines link the different results obtained for the same source 
galaxy when changing the PSF description. The dashed line indicates the 
limit below which galaxies are considered too small to be used in the Eu- 
clid weak lensing analysis. 



Figure [6] also shows that for all galaxies in Table [3] the bias 
is larger for the reference Gaussian PSF (black points) than for the 
more realistic PSF3 (blue diamonds). However, the behaviour of 
the bias as a function of FWHM/FWHMpsf is very similar for 
the two PSF models. In particular, the bias is still a few percent 
when FWHM ~ 1.25 FWHMpsf- Note, however, that in the 
case of a PSF as small as that of Euclid we do not expect galax - 
ies with observed sizes similar to the PSF (e.g. lMassev et alj2012l) . 
For this reason, a characteristic size of about 1.4 FWHMpsf can 
be considered representative for the Euclid galaxy sample. For such 
galaxies we expect a bias of a few times 10 -3 . 



5 CALIBRATION OF THE BIAS 

In order to obtain accurate constraints on cosmological parameters, 
the multiplicative bias fo r the Euclid cosmic shear analysis need s 
to be less than 2 x 10" 3 dMassev et all2012MCropper et alj|2012h . 
The results presented in the previous section suggest, however, that 
the spatial variation of the SED will lead to multiplicative biases in 
the ellipticity that exceed the allowed range in the case of Euclid. 
Hence a way to mitigate the problem of colour gradients is required. 

In principle it should be possible to model the colour gradi- 
ents (and thus determine the bias) using resolved images of galaxies 
taken in different filters. This approach was already suggested by 
IVoigtetalj|2rJl2h . Unfortunately, many factors influence the qual- 
ity of the results. The accuracy with which one can model the bias 
depends on the properties of these images: signal-to-noise ratio, 
resolution, pass-band characteristics, number of filters, properties 
and k nowledge of the PSF. 

IVoigt et all d2012l) explored the possibility of a calibration 
sample, covering the full range of properties of the source galaxies, 
for which the bias can be determined. Following their work, we ex- 



amine this route in more detail, starting with the question what data 
are required to determine the bias observationally. 

5.1 Approximated SED reconstruction 

To quantify the colour gradient bias, resolved observations in at 
least two filters are needed. In this section we explore how well 
one can reconstruct the spatial colour distribution of a galaxy when 
observations in only two bands of that galaxy are available. For 
each of the narrower filters, we define the observed image: 



(x)= / AT l (A)5°(x;A)*P,(x;A)dA 



i = l,2. (22) 



We use the observed image in each filter, J° ba (x), to derive the 
approximated intensity S ,0,approx (x; A) that we need to estimate 
the bias. We will for the moment ignore the fact that a galaxy ob- 
served in the narrower filters has already been convolved by a PSF 
Pi(x; A) (i.e. we take it to be a 5-function). This will be addressed 
in the next section. In addition we make the assumption that for 
each pixel the SED can be interpolated linearly: 



-rO.approx 



(x;A) = a(x)A + 6(x). 



(23) 



The coefficients (a, b) can be determined for each pixel by solving 
a linear system of equations: 



AT, 



AA, 



i(A) [a(x)A + 6(x)]dA 



robs / \ 

li (x) 



1,2. 



(24) 



Once we obtain the approximated intensity 5 , 2 pprox (x; A) we use 
it to estimate the bias induced by the colour gradients in the same 
way we have done previously. 

As shown below, HST observations represent the best avail- 
able sample to study the impact of colour gradients on the shear 
estimation, which explains our choice of filters. We use the spatial 
resolution and transmission 21 (A) for the F606W and F814W fil- 
ters (see left panel of Figure [2). Note that this procedure can be 
applied to other sets of resolved observations taken in two or more 
filters, although the accuracy of the results will vary. The use of 
more bands allows one to estimate the SED more accurately, for 
example by using higher order polynomial interpolations. 

In the presence of noise, as is the case for real data, it will 
be more challenging to reconstruct the spatial colour distribution 
by interpolation. In the case of Euclid we do not expect this to be 
a major issue as the galaxies used as calibration sample will have 
a relatively high signal-to-noise ratio in typical HST observations. 
Alternatively, one might consider ways to reduce the im pact of the 
noise , for example b y using a shapelet decomposition (Refregier 
l2003l : lKuiikenll2006l) . 

Having fixed the properties of the two filters, we can compute 
the difference Am between the bias for the Euclid broad-band ob- 
servations and the bias reconstructed using the linear interpolation 
from F606W and F814W observations. Figure [7] shows the results 
for model galaxy B as a function of redshift when we vary the pa- 
rameters as was done for Figure [5] We find that we are able to 
recover the bias with an accuracy of a few times 10 -4 . 

We note that Am generally shows the same features and that 
its amplitude is proportional to the input bias. This stems directly 
from the fact that the linear interpolation cannot reproduce the 
Irr SED accurately enough (see Figure [T}. The main reason for 
this is associated with the presence of emission lines in the range 
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Figure 7. The difference between the predicted bias for an Euclid-tike broad-band filter and the bias obtained by modelling spatially resolved observations in 
the HSTF606W and F814W filters. The results are shown as a function of redshift for model galaxy B, varying the same parameters as was done in Figure[5] 



400 nm < A < 600 nm, although this is not the only reason. Not 
surprisingly, the linear interpolation cannot reproduce any of the 
SEDs perfectly over the large range in wavelength covered by the 
Euclid pass-band. The linear interpolation fails to capture some of 
the features of the SEDs visible at A < 500 nm and the Balmer 
break at 400 nm. 

Hence the inaccuracy in modelling the SED at each position 
leads to a residual bias which is a strong function of redshift. How- 
ever, we note that Figure[7]exacerbates the problem because the lo- 
cal Irr SED is always the same and the linear interpolation fails co- 
herently. In practice, the linear interpolation will sometimes over- 
estimate and sometimes understimate the bias such that the average 
bias as a function of redshift might still be estimated correctly and 
the residuals uncorrelated. The SED also depends on the age of the 
stellar population, the metallicity, the dust extinction, and the ve- 
locity dispersion; all aspects we have neglected. For instance, the 
bottom-right panel of Figure [7] shows that changing the Irr SED 
used for the disk to any of the Sa-Sd spectra changes the residual 
bias substantially. 

Finally, it is possible to improve on these results by perform- 
ing an SED fit instead of a linear interpolation. This might be feasi- 
ble especially in the case of low redshift galaxies for which multi- 
colour observations are available and can be used to put tight con- 
strains on the local SED. This could significantly improve the ac- 
curacy on the estimated bias at low redshift and help to reduce the 
residual bias at higher redshifts. 



5.2 Effect of the native PSF 

Since the images we use to evaluate the colour gradient-induced 
bias have been convolved with a PSF themselves, the procedure 
used to retrieve the local SED is more complicated in practice. Ac- 
counting for the PSF, Equation d22t changes to 



/ Ar i (A)[a(x)A + 6(x)]*P I (x;A)dA = / l obs (x) , i = l,2. 

J AAj 

(25) 

In Fourier space we obtain a linear system of equations which we 
can solve to obtain a(k) and b(k), the Fourier transformed maps 
of the linear coefficients to approximate the local SED. Solving a 
linear system in Fourier space corresponds to performing a decon- 
volution. This will therefore cause loss of information, as we cannot 
reconstruct scales that are smaller than the PSF. Note that this sys- 
tem of equations can be expanded to more filters, and higher order 
interpolations. The filters might also have different PSFs. 



6 EUCLID BIAS MODELLING WITH HST FILTERS 

In the previous section we have shown that it is possible to quan- 
tify the effect of colour gradients using resolved observations in at 
least two filters. A complication is that one needs to deconvolve the 
galaxies to account for the native PSF of the narrower filters. This 
will reduce the overall accuracy of the measurement of the colour 
gradient bias for two reasons. First of all, solving the system of 
equations {25} for the local SED is equivalent to a deconvolution, 
implying that there is an upper limit to the spatial frequencies we 
can recover. Secondly, the native PSF we need to correct for is not 
perfectly known. 

Since deconvolution algorithms have intrinsic limitations 
which depend on the size of the PSF, a small PSF is always prefer- 
able. Amongst currently available data, HST observations are there- 
fore the most suitable to model colour gradients. As can be seen in 
Figure [8] the characteristic size of the Euclid PSF is twice the HST 
PSF. This stems directly from the fact that the diameter of the Eu- 
clid mirror D — 1.2 m is about half of the size of the HST mirror, 
D = 2.5 m. Despite its small size, there is a loss in accuracy caused 
by the HST PSF. Furthermore, the HSTFSF does vary as a function 
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Figure 8. Top panel: Comparison of the Euclid PSF profile at A = 550 nm 
(solid-dotted line) with the ACS (Wide Field Channel)-rtST PSF profiles at 
the same wavelength. For the HST PSF we show the different profiles as 
a function of position and defocus values expressed in fim. Bottom panel: 
the same as top panel but for A = 800 nm, 



of position and time, which can be modelled with finite accuracy. 
The impact of both complications are quantified below. 

We generate images of the HST PSF at different positions 
and for various fo cus configurations using the Tiny-Tim software 
dKrist et alj|201 ]]). The PSF images P(x; A) have a pixel size of 
0.05 arcsec and are sampled as a function of wavelength with a step 
<5A = 50 nm. When convolving the images of the source I°(k; A) 
with the HST PSF, we approximate the PSF to the closest one in A 
instead of interpolating between the various PSF images. 

The top panel of Figure [8] shows our nominal HST PSF pro- 
file, which is the one in the middle of the first ACS chip, without 
defocus, for A = 550 nm. We compare this profile to those ob- 
tained by changing the position across the camera and the focus. 
In the bottom panel we show the same HST PSF profiles but for 
A = 800 nm. Note that we ignore the camera distortion, as it does 
not alter colour gradients and can be corrected for before modelling 
the bias. 

Both focus and position changes from the nominal configura- 
tion increase the half-light radius of the PSF, but in different ways. 
The change of the focus (upper panel of Figure [8} does not affect 
the core of the PSF but it affects the wings, which become larger 
(effectively the PSF profile is slightly suppressed at small scales 
and boosted at larger scales). Changing the position affects both 
the core and the wings of the PSF. This is the result of a change in 
the diffusion coefficients which increases the overall characteristic 
size of the PSF. In both cases the effect is very small and we do not 
expect the final results to depend significantly on the position nor 
on the focus. 



6.1 PSF deconvolution 

We use the PSFs generated with Tiny-Tim for the nominal focus 
and position to create images for model galaxies B and S. We esti- 
mate the bias without accounting for the PSF and present the results 
in Figure [9] (red lines). Comparison with the true bias (black lines) 
shows that ignoring the PSF leads to an underestimate of the bias 
because the native PSF blurs the colour gradients. As expected, the 
effect is larger when the galaxy is smaller (right panel). 

If we instead use Equation J25b to reconstruct the approxi- 
mated intensity profiles we obtain the blue solid lines in Figure [9] 
Hence, when the convolution by the HST PSF is accounted for, we 
recover the original bias quite well, as shown in the bottom pan- 
els where we plot the difference Am between the original bias and 
the estimated one. For the galaxy B we are able to recover the bias 
within an accuracy of 2 x 10 -4 ; as expected, the performance is 
worse for the S galaxy. This difference includes the error made by 
approximating the local SED with a linear function and it is inter- 
esting to note that the residuals for both galaxies are similar to the 
one estimated in Section [5] where we ignored the PSF of the nar- 
rower band data. This suggests that the linear interpolation of the 
SED is still the main limitation and that HST data are suitable to 
model the bias induced by the presence of colour gradients. As a 
caveat we note that the deconvolution step will be more compli- 
cated in practice due to the presence of noise in real data. 



6.2 Imperfect PSF knowledge 

The HST PSF cannot be modelled perfectly because it varies with 
time. The resulting error in the PSF model will lead to an error in 
the estimate of the colour gradient bias. We estimate the impact 
of this by examining the variation in the bias for a range of HST 
PSFs. Under the assumption that the deconvolution and the SED 
interpolation work perfectly, the loss of accuracy is given by the 
dispersion between the bias estimates for the various PSFs. This 
allows us to estimate an upper limit to the accuracy due to the errors 
in the HST PSF model. 

To do so, we evaluate the bias for slightly different PSFs but 
without deconvolving the PSF. We generate PSFs using Tiny Tim 
for various focus configurations and camera positions and measure 
the resulting colour gradient bias using Equation i24l . The top pan- 
els of Figure [10] show the amplitude of the bias for a few positions 
and defocus values. For reference we also show the true bias (black 
solid line). The difference between the bias for the reference HST 
PSF and the bias obtained when the defocus and the position are 
changed is presented in the bottom panels. Note that the differences 
are always smaller that 10 -4 and always positive. This is because 
the changes in position and focus both increase the PSF size (and 
thus the absolute value of the bias). As expected, the impact of 
variations of the HST PSF are larger for galaxy S. For galaxy B, the 
changes in the bias values are about 10~ 5 , whereas for the galaxy 
S they are about 10 -4 in the worst case. 

Were it possible to perform a perfect deconvolution, ignor- 
ing PSF variations as a function of position and focus would lead 
to an extra error term. This error can be roughly estimated as the 
product of the dispersion from the various focus and position con- 
figurations, multiplied by the ratio between the true and recovered 
bias (black and orange lines, respectively). In practice the PSF as a 
function of positi on and focus can b e determined reasonably well 
(see for example iRhodes et alj|2007l ; [Schrabbac k et alj|2010l) and 
the error is expected to be significantly smaller than what is shown 
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Figure 9. Left figure, top panel: Comparison of the true bias for model galaxy B (black solid line) with the bias using a linear interpolation of the SED using 
the F606W and F814W filters, but ignoring the effect of the HST PSF. The red solid line indicates the resulting bias for the reference HST PSF. The blue solid 
line shows the bias estimate when we do account for the HST PSF. The bottom panel shows the difference between the true bias and its estimate accounting 
for the HST PSF. The right panels show the same as the left panels but for the reference galaxy S. 



in Figure [Tol Hence, the limited accuracy of the model for the HST 
PSF is not an important source of error. 



7 CALIBRATION OF THE BIAS 



IVoigt et alJ d2012n proposed to use HST observations of galaxies 
to determine the mean bias as a function of galaxy properties. The 
precision with which this can be done depends on the intrinsic vari- 
ation in galaxy properties, which in turn drives the size of the sam- 
ple of galaxies that is needed. In this section we examine whether 
the HST archive contains a sufficient number of resolved galaxies 
observed in at least two bands to model the bias with the precision 
required for Euclid's science objectives. 

Throughout this paper we have made realistic but rather con- 
servative assumptions in order to estimate the bias using simulated 
bulge plus disk galaxies that by design showed significant colour 
gradients. The actual amplitude of the bias, however, needs to be 
derived using actual observations. A preliminary an alysis based 
on 12 ,000 galaxies in the Extended Groth Strip (EGS. lDavis et alJ 
1 20071) suggests an average value of (m) ~ 3 x 10~ 3 in the worst 
cases, with an uncertainty <7{ m ) of a few times 10 -4 computed in 
redshift bins each containing about 500 galaxies (Huang et al. in 
prep.). 

The average bias is larger than the 3 x 10~ 4 measured by 
IVoigt et alJ J20E% using their full dataset comprising ~ 700 galax- 
ies. A source for this difference is our use in this paper, and in 
the preliminary analysis of EGS data, of a PSF that has a larger 
size and a stronger A dependence than the one used bv lVoigt et alJ 



fe012h which is closer to Euclid's one. Additional differences might 
arise from the fact that we use different shear estimates. Further- 
more, as pointed out bv lVoigt et al]fe012h . their selection of galax- 
ies may not be very representative. For instance, the EGS dataset 
may contai n a fraction of galax ies with red bulges larger than the 
on es in the ISimard et alJ J2002h catalogue; as shown in Figure 6 
of lVoigtet ail d2012h . se lecting galaxies w ith a redder bulge en- 
hances the bias to 10 " 3 . IVoigtetal1 ( l2012h estimate a dispersion 
in the bias of ~ 7 x 10~ 3 for the whole sample, which suggests 
that about a thousand galaxies per redshift bin are needed to ob- 
tain <77 m \ — 2 x 10~ 4 . The error on the average we obtain from 
the preliminary EGS analysis agrees with these values. As shown 
in Table [4] the HST archive contains enough galaxies to reach this 
precision. 



7.1 Residual bias correlations 

The bias depends on the intrinsic properties of the galaxies, which 
in turn depend on environment. The colour gradient bias will there- 
fore vary spatially. A simple correction using the average bias for 
each tomographic bin may therefore not be sufficient, and still lead 
to a spurious signal in the two-point cosmic shear statistics. 

We assess the residual bias on the correlation function as fol- 



3 IVoiet et al l JioTH) quote a number of galaxies that is ten times larger to 
account for possible limitations in the accuracy of their analysis, which is 
based on a small sample of galaxies with a bulge and disk decomposition. 
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Figure 10. Left figure, top panel: Comparison of the true bias for model galaxy B (black solid line) with the bias using a linear interpolation of the SED using 
the F606W and F814W filters, but ignoring the convolution by the HST PSF. The orange solid line indicates the resulting bias for the reference HSTPSF. The 
other lines indicate how the bias changes as a function of the position and defocus value expressed in fj,m. The bottom panel shows the difference between the 
value of the multiplicative bias obtained for a galaxy observed using the reference HST PSF, and the bias when the position of the galaxy is changed and/or a 
defocus is present. The right panels show the same results, but for the S galaxy. 



lows. We compute the two-point ellipticity correlation function 
£+ (8) selecting sources belonging to the same redshift bin: 

= ( 7 t(0%(0 + e')) + ( 7r (0 , )>(0 + e')) , (26) 

where the ensemble average is meant as an average over all pairs 
with distance 9 and 74 and j r indicate the tangential and the 
45-degree rotated components of the estimated shear projected 
along the line connecting the pair of galaxies (see for example 
iBartelmann & Schneiderll200ll) . We define A£(0), the relative bias 
in the correlation function: 



(27) 



Based on the breakdown presented in ICropper et all J2012I) 
we assume that Euclid's science objectives can be achieved if the 
value |A£((?)| is smaller than 5 x 1CP 4 for all angular scales 
6 > 6 arcmin. This limit allows us to estimate the required min- 
imum size of the calibration sample that is used to determine the 
average colour gradient bias as a function of galaxy properties and 
redshift. 

Since we lack observational results we use simulations to con- 
servatively estimate the impact of the spatial variation in the colour 
gradient bias. In the following we use model galaxy S, which has 
a bias of a few times 10 -3 and an observed FWHM of about 
1.4 FWHMpsp, typical of a galaxy in the Euclid survey. 

To simulate the properties of galaxies in the Euclid survey and 



the calibration sample we use ra y tracing simulations produced us- 
ing th e Millennium Simulation (Hilbert et al. 2009; Sprin geTet alj 
2005). The ray tracing simulations comprise 32 lines of sight, each 
covering 16 deg 2 . For each galaxy we know the shear, the redshift 
and magnitudes in various SDSS bands. To simulate the depth of 
the Euclid catalogue we add the i, r and a third of the 2 fluxes of the 
SDSS bands and select galaxies brighter than m r i z = 24.5. The re- 
sulting density is 32 gal/arcmin 2 , in agreement with expectations. 
We use the SDSS r — i values to define the colour of each galaxy. 



7.2 Size of the calibration sample 

We searched the HST archive for ACS observations in at least 
two (suitable) filters. The results are listed in Table [4] where we 
note that we also include scheduled observations. The sample con- 
tains approxi mately 1536 arcmin 2 of F606W and F814W data in- 
cluding EGS jDavis et al.1 12007|) and the o verlap of CANDLES 
dGroginetal]|201ll) with either CO SMOS dScoville et ail |2007|) 
and GOODS (Giavalisco et al.ll2004l) . We have also included the 
area of GEMS dRix et alj|2004l) data that does not overlap with ei- 
ther GOODS or CANDLES. For this part of the GEMS data we 
have F850LP observations which are shallower and slightly redder 
than the F814W data; these can be also used to calibrate colour 
gradients, but the lower signal-to-noise is expected to reduce the 
accuracy. This results in a total area of 2056 arcmin 2 containing 
about 62, 000 galaxies for which we can determine the colour gra- 
dient bias observationally. 
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Figure 11. Left panel: Distribution of 'observed' SDSS r — i colours of the simulated galaxies in three redshift bins. For each redshift bin, the blue (red) 
solid line represents the absolute value of multiplicative bias predicted by our model (see text) as a function of observed colour for the sources with the lowest 
(highest) redshift included in that bin. Right panel: Average bias (m) and its error a({m)) using the calibration sample of ~ 62, 000 galaxies to which we 
assign a bias as a function of colour and redshift as described in Section|7] 



Name Area (arcmin 2 ) Number of galaxies 



EGS 650 18000 

CANDELS/UDS 198 6000 

CANDELS/COSMOS 198* 6000 

CANDELS/GOODS-CDFS 300* 9000 

CANDELS/GOODS-NORTH 190* 5500 



GEMS/CANDLES+GOODS 520 (F850LP) 15500 



crease in its flux. The galaxy S has a bulge size of th = 0.09 arcsec 
which contains 25% of the light; we set ru of the bulge to that of 
the disk (r^ = 0.59 arcsec) when the flux of the bulge is 100% 
(note that in this case we have an elliptical galaxy with no colour 
gradients) and create all other cases using a linear relation between 
flux and size. 

Once we have created this set of model galaxies, we derive 
their colours by integrating the flux over the SDSS i and r fil- 
ters. The resulting colours are generally too blue compared to the 
galaxies in the Millennium Simulation. These differences are likely 
caused by differences in the adopted SEDs. The change in colour 
with redshift for our model early-type galaxies is the same as that 
of the red Millennium galaxies. Thus, we correct for this offset in 
colour by shifting them slightly to match the results from the Mil- 
lennium Simulation. 

The left panel of Figure QTJ shows the colour distribution of 
the simulated Euclid-like galaxies for three redshift bins. For each 
redshift bin we also show the bias as a function of the observed 
colour. The blue (red) lines indicate the bias for a galaxy at the 
lower (upper) limit of the respective redshift bin. In all cases the 
bias reaches a maximum value when the bulge contains about 20- 
30% of the light. 

The bias m assigned to each galaxy in the Millennium Simu- 
lation is computed by interpolating the value of the multiplicative 
bias m as a function of the r — i colour at the redshift of the galaxy. 
The average bias as a function of redshift is presented in the right 
panel of Figure QTJ Figure QTJ shows that —4 x 10~ 3 < m < 
—2 x 10 -3 , similar to what is found from the preliminary EGS 
analysis. The dispersion of the simulated bias, however, is about 
1 x 10~ 3 , smaller than that found from the EGS analysis. To ensure 
that our estimates remain conservative, we increase the variance of 
m accordingly: we include an additional spread in the bias sampled 
by a Gaussian of width a — 6 x 10 -3 . 

This results in a simulated data set where the bias depends 
on redshift and colour with an extra variance to account for varia- 
tions as a function of other parameters (such as ellipticity and size). 
These estimates were used to compute the errors on the mean bias 
as a function of redshift presented in Figure QTJ where we assume 



total 2056 62000 

Table 4. Size of the HST data sample observed in both F606W and 
F814W/F850LP bands. The entries marked with * are not yet (fully) avail- 
able, but will be soon. All observations quoted in the table are deeper than 
Euclid. The number of galaxies has been computed assuming F814W < 
24.5 which matches the number density of 30gal/arcmin 2 expected for 
Euclid. Note that the F850LP observations are shallower but they can still 
be used with in combination with the F606W data to obtain an estimate for 
the local SED. 

The data listed in TableUJwere obtained as part of several sur- 
veys. Including single ACS pointings as well would increase the 
area by another ~ 1000 arcmin 2 . Finally, d edicated deep Eu clid 
observations of the area covered by STAGES dGrav et alj 2009) and 
COSMOS would increase the sample by more than a factor of three. 
The full benefit of the latter observations requires more study be- 
cause the Euclid PSF size is larger and pass-band broader. In the 
following we therefore conservatively assume that we can measure 
the colour gradient bias using HST observations of 62, 000 galax- 
ies. 

7.3 Bias model 

To each galaxy in the simulation we assign a bias m, which is a 
function of colour and redshift. Including the colour dependence is 
important to capture the fact that galaxy colour depends on envi- 
ronment. The value for the bias is obtained by taking the S galaxy 
and varying the fraction of light in the bulge as we did in the top- 
left panel of Figure [5] Note that increasing the fraction of light in 
the bulge without changing its size results in redder galaxies with 
an unrealistically small but very bright bulge. To avoid this, we in- 
crease the size of the bulge accordingly, thus accounting for the in- 
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Figure 12. Left panel: The spurious contribution to the ellipticity correlation function due to colour gradient bias £ c ($) (red solid line) compared to the cosmic 
shear signal £+(8) for a tomographic bin with 0.2 < z < 0.3 (black solid line). The magenta dashed line shows the angular scale corresponding to the Euclid 
field-of-view. The green triangles indicate the spurious signal if the bias correction is a function of redshift only. If the correction is a function of both colour 
and redshift the blue triangles are obtained. Right panel: The absolute value of the relative change in the ellipticity correlation function if £ c = 0. The colours 
and symbols are the same as for the left panel. The dotted black line indicates the maximum allowed value for [ A£| for Euclid to meet the requirements. The 
eiTors represent the standard deviation from 32 independent realisations of the calibration sample (see text). 



a calibration sample of 62, 000 galaxies (40, 000 of which have 
z < 1.2). 



7.4 Effect on the ellipticity correlation function 

So far we have limited the discussion to the impact of colour gra- 
dients on the multiplicative bias m. As the bias is related to the 
correction for the PSF, we expect the additive bias c to be affected 
as well. In this more general case, the shear estimate for the i- 
component is given by: 



7< 



(1 + 771)7; + c - 



Inserting this expression into Equation j26\ we find that the mea- 
sured correlation function can be expressed as: 



■■Z+{0)[l + 2(m)+U(0)] +&(*): 



(28) 



where £ m {0) and £. c {9) quantify the correlations of m and c, re- 
spectively, as a function of angular separation. The relative bias in 
the correlation function is thus given by: 



A£(0) = [2 (m)+£ m (0)] + 



gc(g) 



(29) 



lMassevetaf](l2012h derived expressions for the various sources of 
multiplicative and additive bias and showed that the contributions 
to m and c caused by colour gradients are related through: 



epsp.i 

1 Tp — 

1 7 i epsF 



(30) 



wh ere epsF.; is the i- component for the PSF ellipticity. Follow- 
ing [Massey^etal] J2012h . we assume -P epSF = 1. Hence we have 



£c(0) = U((?)£+,psf(0)/P 7 2 , where £+,psf(0) is the elliptic- 
ity correlation function of the PSF anisotropy. Because our defi- 
nition of ellipticity b y Equation U5i differs from the one used by 
iMassev et al.l d2012l) . we have P 1 ~ 1 — e 2 ~ 0.86 whereas we 
adopt epsF = 7% for the PSF ellipticitjQ, which is an extremely 
pessimistic scenario. 

Although c is smaller than m, Equation l |29t shows that the 
additive bias can become dominant on large scales if the PSF el- 
lipticity correlation function does not vanish sufficiently quickly. 
In our case the PSF orientation is constant across the field, which 
implies that for any given value of c there will be a scale for which 
(,c{d) > £+(0), unless £ m (6*) vanishes faster than This scale 

is smaller for low-redshift sources for which the cosmic shear sig- 
nal is smaller. 

This large-scale signal provides a way to estimate the ampli- 
tude of the colour-gradient induced bias directly from the Euclid 
dataset by measuring the correlation between the ellipticity of stars 
and the PSF-corrected ellipticity of galaxies (i.e. the shear). This 
correlation is commonly used as a di agnostic to ensure th at the 
PSF anisotropy is well corrected (e.g. iHevmans et all l2012t> . For 
this reason one would rather not use this residual signal to correct 
for colour gradients. However, in the case of colour gradients we 
do not need to correlate galaxies and stars in the same field: the 
colour gradient bias will be proportional to the PSF anisotropy, so 
that residuals of independent fields will be correlated. Note that in 
the absence of colour gradients these correlations are supposed to 
vanish if the correction for PSF anisotropy is accurate on average. 

The calibration sample is used to obtain an estimate for (m), 
the average value of the bias for a particular selection of sources, 
which is used to correct the shear signal for these galaxies. Hence 

4 This corresponds to the maximum value for the polarisation of 15% 
allowed for the Euclid PSF lLaureiis et al.l201 ltlCropper et alj20l"3) . 
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redshift only redshift and colour 

0.2 < z < 0.3 0.5 < z < 0.6 0.8 < z < 0.9 0.2 < z < 0.3 0.5 < z < 0.6 0.8 < z < 0.9 



8 < 240 arcmin 34.5 ± 14.6 1.0 ±11.8 2.1 ±5.4 5.7 ±4.4 -5.1 ±11.3 -2.2 ±2.2 

0>6arcmin 31.3 ±16.1 -1.4 ±13.4 1.3 ± 6.5 7.4 ±4.8 -6.3 ±14.2 -2.8 ±2.6 

6 < 9 < 50 arcmin 28.0 ± 8.9 -1.7 ±1.5 1.2 ± 0.4 6.1 ± 1.6 -3.6 ±3.3 -1.3 ±0.3 

Table 5. Average relative bias A£(#) in units of 10~ 4 for the three redshift bins presented in Figure [T3l Columns 2 — 4 list results if the bias is modelled as a 
function of redshift only; Columns 5 — 7 list result when the bias is modelled as a function of both redshift and colour. 
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Figure 13. Absolute value of the bias A§(0) as a function of scale for three redshift bins. The dotted line shows the maximal residual bias allowed to meet 
Euclid's science requirements. For reference, the vertical dashed line indicates the size of the Euclid field-of-view. The red line shows the relative bias in 
the ellipticity correlation function when colour gradients are ignored. In this case the multiplicative bias lowers the signal (since m < 0), but the additive 
bias leads to an increase in signal. The green triangles show the relative bias when the correction is a function of redshift only. The results improve when the 
correction is a function of colour and redshift (blue squares). The errors represent the standard deviation from 32 independent realisations of the calibration 
sample. The mean values of the residual bias and their errors are listed in Table[5] 



the relevant quantity becomes 5m, the residual bias for each galaxy 
after correction. After removing the average bias, we can rewrite 
Equation j29\ to obtain: 

Ag(fl) = ft m (0) x (l + ± ) , (3D 

where (,sm{9) is the correlation function of the residuals 
(5m(8)5m(8 )). As for the uncorrected case, the first term is re- 
lated to the multiplicative bias, while the second term corresponds 
to the additive bias. 

7.5 Bias estimation 

We assume the calibration sample to be a contiguous patch with the 
same area as the total area of the datsets presented in Table|4] This 
increases the sampling variance compared to what one would ob- 
tain from several independent fields. Hence, our estimates are con- 
servative. The model for the bias from this calibration subsample is 
then used to correct the bias for the whole survey. To estimate the 
sampling variance we generate 32 independent realisations where 
each time a different location for the calibration sample is selected. 

We create simulated shear catalogues for three different sce- 
narios. In the first scenario we do not correct for the bias, whereas 



in the second case we use the calibration sample to derive an aver- 
age bias value as a function of the redshift. We sample the bias in 
steps of Az = 0.1. In the third case we use the calibration sample 
to determine the average bias as a function of redshift and observed 
colour. We sample the bias in steps of Az = 0.1 and 15 bins of 
colour of width 0.10 covering the range between [0.05, 1.45]. 

We start by exploring the m and c contributions to the bias 
in the correlation function separately. The left panel of Figure [T2l 
shows the amplitude of £ c (#) for the tomographic bin 0.2 < z < 
0.3, which can be compared to the cosmic shear signal (6) (black 
line). Since the PSF orientation is constant, f c (0) is always positive, 
its asymptotic value is (m) 2 ef> SF and there is a scale for which the 
value of £ c becomes larger than the cosmological signal. This hap- 
pens at scales larger than 200 arcmin even for the smallest redshift 
bin. In reality £ c (#) is expected to be a factor ~ 10 smaller as 
the Euclid PSF anisotropy in the currently agreed optical design (c. 
201 1-2012) is supposed to be about 7% only in the corners while it 
is about 2 — 3% in most of the field of view. Additionally, the PSF 
pattern varies across the field, so that £psf(#) decreases quickly 
and it is expected to be very small on scales that are larger than the 
field of view (indicated in the left panel of Figure[l2]and in Figure 
I13lbv the magenta dashed lines). 

As indicated by the green triangles in Figure [T2] the value for 
£ c (#) is reduced by two orders of magnitude at large angular scales 
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when we use the average bias as a function of redshift to correct the 
colour gradient bias. The improvement is smaller on small scales, 
where the clustering of galaxies is important. This is remedied by 
modelling the bias as a function of colour and redshift (blue trian- 
gles). The error bars shown in Fi sure [T2l have been computed from 
the 32 realisations and include sampling variance. 

The right panel of Figure[T2]shows the value of A£(0) when 
£ c (#) = 0. As expected, when no correction is applied, the overall 
bias |A£(#)| is about 2 (m) |. Since m < 0, the signal is sup- 
pressed (i.e A£(#) < 0). Once again, correcting the average bias 
as a function of redshift improves the results. Including the colour 
dependence improves the results even further to a level below the 
requirement of 5 x 10 -4 (indicated by the dotted line). 

Figure Q~3] shows \A£(9)\, for three redshift bins. If the colour 
gradient bias is ignored, | A^^) | has a minimum due to the compe- 
tition between the multiplicative and additive terms, which have an 
opposite effect on the observed ellipticity correlation function. The 
location of the minimum indicates the scale where the additive bias 
starts to dominate. Similarly to what we found in Figure[l2] correct- 
ing the measured shear by the average bias as a function of redshift 
reduces the amplitude of the residual bias on the correlation func- 
tion. Taking into account the dependence upon the observed colour 
reduces this residual further, in particular for small angular scales. 
For reference we indicate with a dotted line the maximum value for 
A£ (9) | allocated to achieve Euclid's science objectives. 

The measurements at different angular scales in Figure[73lare 
correlated, which complicates a simple interpretation of the sig- 
nificance of the results. We therefore computed the average resid- 
ual bias, accounting for the correlation between the angular scales 
using the covariance matrix computed from the simulations. The 
results are reported in Table [5] for the three redshift bins after cor- 
recting the bias either as a function of redshift or as a function of 
both redshift and colour. We derive the average residual bias for 
three different ranges: the full range of scales, i.e. 9 < 240 arcmin; 
eliminating scales which will not be included in the cosmological 
analysis, i.e. 8 > 6 arcmin; eliminating scales larger than the Eu- 
clid field-of-view, i.e. 6 < 9 < 50 arcmin. In all cases we find that 
modelling the bias as a function of colour and redshift produces 
better results than modelling the bias as a function of redshift only. 
With this correction we obtain residual biases that are consistent 
(within la) with the maximum allowed value (i.e. 5 x 10 -4 ). A 
more realistic model for the PSF anisotropy would result in a sig- 
nificantly smaller residual bias, especially on large scales. 

Overall these results demonstrate that the HST archive con- 
tains a sufficient number of galaxies observed in at least two suit- 
able bands to calibrate the colour gradient-induced bias to the pre- 
cision required for Euclid. To simulate the spatial variation of the 
bias using the Millennium Simulation, we assumed that the colour 
is the main parameter to consider. We thus ignored the fact that the 
bias also depends on other parameters such as its ellipticity, Sersfc 
index, etc., as is evident from Figure|5] This has the effect of reduc- 
ing the simulated scatter in the bias for a given colour and redshift, 
but we accounted for this by including an intrinsic scatter in order 
to match results from a preliminary analysis of EGS data (Huang 
et al., in prep.). We note that we assumed a rather pessimistic PSF 
anisotropy and PSF size and therefore are confident that our esti- 
mates are in fact conservative. On real data, we can can test further 
bias dependencies and model it as a function of an optimal set of 
parameters. The efficiency of the correction can be tested and even- 
tually improved by dividing the Euclid dataset in subsamples and 
comparing the differences in shear distributions measured on these 
subsamples. 



8 CONCLUSIONS 

Euclid wi ll measure the cosm ic shear signal with unprecedented 
precision jLaureiis et al.ll2oTll) . To do so it will measure galaxy 
shapes from images observed in a rather broad filter. In combina- 
tion with the wavelength dependence of the PSF this leads to biases 
in the shape measurements because the colours of galaxies vary 
spatially. To ensure that the constraints on cosmological parameters 
are not compromised it is important to quantify this 'colour gradi- 
ent' bias and examine possible approaches to mitigate the problem. 

In this paper we show how the spatial variation of the colour 
of a galaxy generally leads to a bias in the measurement of its 
shape. In the case of unweighted moments the PSF can be removed 
without introducing any bias. In practice, in order to suppress the 
noise, shear measurement methods use a weight function that as- 
signs more weight to the inner regions. Consequently, the bias de- 
pends on the choice of the weight function and is therefore gener- 
ally method-dependent. We derive results for moment-based meth- 
ods, but note that our findings are also applicable to methods that 
fit galaxy models to the data. 

We showed analytically that the amplitude of the bias scales 
proportional to the square of the width of the filter used for the mea- 
surement. We studied in detail how the bias depends on the char- 
acteristics of the galaxies and the PSF using model galaxies. The 
amplitude of the bias depends strongly on the ratio of the observed 
FWHM of the galaxy to the FWHM of the PSF, with a larger PSF 
leading to a larger bias. The bias can be as large as a few percent if 
the observed galaxy is about the size of the PSF, while it decreases 
rapidly to only a few times 10~ 4 when the galaxy size is 1.6 times 
the PSF size. Given the small size of the Euclid PSF and the mag- 
nitude limit of 24.5, the observed galaxies will generally have an 
observed FWHM which is larger than the PSF. The bias depends 
on the bulge and disk sizes, on the ellipticity, but it depends most 
strongly on the colour of the galaxy. 

It is possible to determine the bias observationally using spa- 
tially resolved data observed in at least two filters. We focus on the 
particular case of HST observations in F606W and F814W, which 
are best suited for this purpose. We quantify the limitations in mod- 
elling the local SED using only two filters taking into account the 
HST PSF. We find that the linear interpolation we use to approx- 
imate the local SED is the main source of error. This could be 
improved in principle by performing a local SED fit. This might 
be possible especially for low redshift galaxies where we have re- 
solved observations in various colours. Errors in the model for the 
HST PSF have a very minor impact. Our results therefore indicate 
that we can model the bias due to colour gradients using HST obser- 
vations with an accuracy of a few times 10 -4 for a typical galaxy 
used in the Euclid weak lensing analysis. 

To determine the average bias with sufficient precisi on re- 
quires a representative sample of galaxies. IVoigt et all d2012h sug- 
gested that a sample of about thousand galaxies per redshift bin 
would be sufficient to determine the bias with a precision of 
<T( m ) = 2 x 10 -4 . This has been confirmed by a preliminary anal- 
ysis using EGS data (Huang et al. in prep.). The amount of archival 
HST data observed in at least two bands is sufficient for this. 

A complication arises from the fact that the bias depends 
strongly on the colour. As a result the clustering of galaxies will 
lead to spatial variations in the bias. In this case an average cor- 
rection for each redshift bin may still lead to biases in the two- 
point ellipticity correlation function. To examine the impact of this 
second order effect we used a ray-tracing catalogue based on the 
Millennium Simulation. We find that correcting the bias as a func- 
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tion of both redshift and colour using a contiguous patch containing 
about 40,000 galaxies is sufficient to correct the bias to the accu- 
racy required for the Euclid two-point shear tomography. This is a 
conservative estimate of the number of galaxies in the HST archive 
at the time of Euclid's launch. We therefore conclude that the pres- 
ence of colour gradients in galaxies is not a limiting factor for the 
Euclid cosmic shear analysis. 
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